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After removing the double-well potential trapping two initially independent Bose condensates, 
, the density expectation value is calculated when both the exchange symmetry of identical bosons 

and interatomic interaction are considered. The density expectation value and evolution equations 
' are obtained based on both the first-quantization and second-quantization methods. When the 

interatomic interaction is considered carefully, after the overlapping of two initially independent 
condensates, it is shown that there is a nonzero interference term in the density expectation value. 
It is found that the calculated density expectation value with this model agrees with the interference 
pattern observed in the experiment by Andrews et al (Science 275, 637 (1997)). The nonzero 
. interference term in the density expectation value physically arises from the exchange symmetry 

(£) • of identical bosons and interatomic interaction which make two initially independent condensates 

become coherent after the overlapping. For two initially independent condensates, our researches 
I , show that there is an interaction-induced coherence process. 

-«-> , 

^ . I. INTRODUCTION 

0: 

I \ The coherence property plays an essential role in the wave nature of Bose-Einstein condensates (BECs) , and it has 
been investigated intensively [l[ after the experimental realizations of BECs in dilute gases 0, H, 0|- For a single 
BEC which is a perfect quantum fluid at zero temperature, the existence of a macroscopic wave function (or order 
parameter) means that spatial coherence is an intriguing property of the condensate. In the dilute Bose-condensed 
gases, the spatial coherence can be investigated directly by interfering two Bose condensates. For two coherently 
separated BECs, because the relative phase of two sub-condensates is locked, it is not surprising that there is a clear 
interference pattern when the two sub-condensates are allowed to overlap. In the celebrated experiment by Andrews 
et al however, high-contrast fringes were observed even for two completely independent condensates at an initial 
time. In this experiment, to prepare two initially independent condensates, the dilute Bose gases were evaporatively 
cooled in a double-well potential created by splitting a magnetic trap in half with a far-off blue-detuned laser beam 
{Nl . (See also @). In particular, the height of the external potential due to the laser beam is much larger than the 
t-H ' chemical potential of two separated condensates, which means that the tunneling current can be safely omitted. After 
switching off the double- well potential, the two initially independent condensates overlapped and high-contrast fringes 
were observed in Ref. Q . This experimental result shows clearly that there is a spatial coherence property after the 
overlapping between two independent interacting condensates. 

Although high-contrast fringes were observed for two initially independent condensates, in many literature (see for 
example [3, Q and references therein) , it is shown that there is no interference term in the density expectation value 
for two initially independent condensates. To solve the contradiction between this result and the observed high-contrast 
fringes [f| for two initially independent condensates, the observed high-contrast fringes were interpreted with the aid 
of the high-order correlation function P (r, r', t) (which is an oscillation function of r — r') and quantum measurement 
process in the present popular viewpoint. Several theories have been proposed to interpret the observed high-contrast 
fringes for two initially independent condensates such as the stochastic simulations of the photon detection for atoms 
[To| , the expansion of Fock state by the linear superposition of coherent states [Tl| , and the continuous measurement 
theory [12]. 

In the present work, we calculate the density expectation value for two initially independent condensates by including 
carefully the interatomic interaction. Quite different from the simple derivation in other theories, it is found that for 
the case of two initially independent condensates, upon expansion, there is a nonzero interference term in the density 
expectation value when the interatomic interaction and exchange symmetry of identical bosons are both taken into 
account. To exclude any possible error in the cumbersome derivations, we calculate the density expectation value 
and evolution equations by both the first-quantization and second-quantization methods, and the same results are 



o 
o 



c3 

S 

i 



X 



* xionghongwei@wipm.ac.cn 



2 



obtained. After removing the double-well potential trapping the two initially independent condensates, we give the 
theoretical result of the density expectation value which agrees with the interference pattern observed in [f| . 

To show the interaction-induced coherence process between two initially independent Bose condensates when both 
the interatomic interaction and exchange symmetry of identical bosons are considered carefully the paper is organized 
as follows. In Sec. II, we give the general expression of the density expectation value for two initially independent 
condensates based on the many-body wave function of the whole system. It is shown that there is a nonzero interference 
term when the wave functions of two initially independent condensates are no more orthogonal after their overlapping. 
In Sec. Ill, we give the general expression of the overall energy of the whole system, and the evolution equations are 
given by the action principle. In Sec. IV, we give a brief proof why the wave functions of two initially independent 
condensates become non-orthogonal after their overlapping when the interatomic interaction is considered. In this 
section, an effective order parameter is introduced to show further the interaction-induced coherence formation process 
for two initially independent condensates. In Sec. V, we give the numerical results for the evolution of the density 
expectation value according to the experimental parameters in Ref. Q . It is shown clearly that our theoretical results 
of the density expectation value agree with the experimental results of the interference patterns. In Sec. VI, the 
evolution of the density expectation value is given for different coupling constants. It is shown that increasing the 
coupling constant has the effect of enhancing the interference effect for two independent condensates. To show further 
the physical mechanism of the interaction-induced coherence process, in Sec. VII, we prove that the theoretical results 
based on the second-quantization method is the same as the results derived from the many-body wave function in 
the previous sections. In Sec. VIII, the role of quantum fluctations and the orthogonality of the many-body wave 
function are discussed. Finally, we give a brief summary and discussion in the last section. 



II. THE EXPRESSION OF THE DENSITY EXPECTATION VALUE FOR TWO INITIALLY 

INDEPENDENT BECS 



First, we give a brief introduction to the scheme of observing the interference effect of two separated condensates. 
As shown in Fig. 1, the double- well potential can be created by superposing a far-off blue-detuned laser beam (which 
generating a repulsive optical dipole force for atoms) upon the magnetic trap. The combined double-well potential is 
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(1) 



where the first term is the trapping potential due to the magnetic trap, while the second term represents the potential 
due to the laser beam. 

For the Bose gases confined in this double-well potential, there are two quite different cases, (i) If the height 
Uq of the external potential due to the laser beam is smaller than the chemical potential of the system, the two 
sub-condensates can be regarded to be coherently separated after the evaporative cooling. In this situation, the 
two coherently-separated sub-condensates can be regarded as a single condensate. The relative phase of two sub- 
condensates is thus locked, (ii) If the height Uq of the laser beam is much larger than the chemical potential of the 
system so that the tunneling effect can be omitted safely, the two condensates are completely independent after the 
evaporative cooling. 

For two coherently separated sub-condensates, every atom is described by the following normalization wave function 



(r,*) = 



/VN, 



(2) 



where <p c i (r, t) and <j) C 2 (r, t) are the normalization wave functions accounting for the two sub-condensates. N — 
N% + N2 is the total number of particles. Before removing the double- well potential, the average particle numbers 
in each condensate are N\ and N2 , respectively. After removing the double-well potential, the evolution of the wave 
function <j> c (r,t) can be obtained based on the Gross-Pitaevskii (GP) equation [H, EH, H3I . The density expectation 
value is then 



n c (r,t) = N\tf> c (r,t)\" 



N x \4> el (r, t)\ 2 + x Re [<j>* cl (r, t) c2 (r, t)]+N 2 |0 c2 (r, t)\' 



(3) 



The second term in the above equation accounts for the interference effect when there is an overlapping between two 
sub-condensates upon expansion. For <j) c (r, i) satisfying the GP equation, the nonlinear effects in the interference 
pattern of two coherently separated sub-condensates were investigated in Refs [HI, H3] • 

As shown in Fig. 1(b), if the intensity of the blue-detuned laser beam is sufficiently high so that the tunneling effect 
can be omitted, the two condensates can be regarded to be completely independent. In this situation, the number of 
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FIG. 1: The scheme of observing the interference effect of two separated condensates. Fig. 1(a) shows the case of two coherently 
separated sub-condensates at an initial time, while Fig. 1(b) shows two initially independent condensates. 



particles Ni and N 2 in each of the two condensates are fixed. We assume that there are Ni particles described by the 
wave function <j)\ in the left well, while there are N2 particles described by the wave function <f>2 in the right well. To 
investigate clearly the role of the exchange symmetry of identical particles and interatomic interaction, we calculate 
in the following the density expectation value (r, t) directly from the many-body wave function. In Sec. VII, we 
will also calculate the density expectation value based on the second-quantization method, and the same result is 
obtained. 

We will prove in Sec. IV that in the presence of the interatomic interaction, tfii (r, t) and <j>2 (r, t) will become non- 
orthogonal after the overlapping between two initially independent condensates. Thus we consider here the following 
general case from the beginning 

C(t) = J 4> 1 (r,t)<p* 2 (r,t)dV=\((t)\e i ^. (4) 

After removing the double-well potential, the two initially independent condensates will overlap, and thus one should 
consider the indistinguishability of identical particles. When the exchange symmetry of identical bosons is considered, 
the many-body wave function is 



P 
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<f>2 (rjVi + i , t)---(t> 2 (y Ni+N2 , t)] 



(5) 



where P denotes the N\/(Ni\N 2 l) permutations for the particles in different single-particle state <fii or <p2- A n is 
a normalization factor to assure J \^n 1 n 2 ( r ij r 2, ■ ■ • ,^N,t)\ dridr2 ■ ■ -cLym = 1. A n is determined by the following 
equation: 



A n 



i=0 



(Nt - i)l (N 2 - i)\ 



1/2 



(6) 



In this paper, to give a concise expression for various coefficients such as A n , we have introduced the rule 0° = 1. When 
the non-orthogonal property between <f>\ and <j) 2 is considered, one should note that the normalization constant A n is 
relevant to the parameter £. From the form of the many-body wave function (|5|), the quantum depletion originating 
from interparticle interaction is omitted. Thus, this form of quantum state is valid when a s /l << 1 with a s and I 
being respectively the scattering length and mean distance between particles. The role of quantum depletion will be 
discussed in Sec. VIII. 

From the above many-body wave function, after straightforward derivations, the exact expression of the density 
expectation value takes the following form: 



n d (r, t) = N J ^* NlN2 (r, r 2) • • • , r N ,t) $ NlN2 (r, r 2 , • • • , r N , t) d 3 r 2 ■ ■ ■ d 3 r N 
= a d |0! (r, t)\ 2 + 2b d x Re [e^ ( r , t ) cp 2 (r, t)] + c d \cf> 2 (r, t)\ 2 , 

where the coefficients arc 



min(AT 1 — 1,N 2 ) 
ad = a d (i) , 

i=0 

min(iVi— ljiVa-l) 

bd = X] h d{i), 

i=0 

min(N 1 ,N 2 -l) 

Cd = ^ Cd ■ 

i=0 



In the above summations, 



a d (i) 
b d (<) 
c d (i) 



A^x!iV 2 !|C(t)| 
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i\i\{N 1 -i-l)\{N 2 -i)V' 



A 2 n Ni\N 2 \\<;{t)\ 



2i+l 



i\ {i + 1)! (iVi —i — 1)! (N 2 -i - 1)1 
i\i\(Ni -i)\(N 2 -i- 1)!' 



(7) 

(8) 
(9) 
(10) 

(11) 
(12) 

(13) 



For two independent ideal condensates, before the overlapping between the two condensates, we have £ (t = 0) = 0. 
Based on the Schrodinger equation, it is easy to verify that after the double-well potential separating the condensates 
is removed, we have £ (t) = at any further time. Thus b d = 0, and the density expectation value is given by 



n d (r , t ) = N x 1 0x (r , t ) | 2 + N 2 \ cf> 2 (r , t ) | 2 . 



(14) 



In this situation, the interference term is zero in the density expectation value. 

In the presence of the interatomic interaction, £ (t) can be a nonzero value after the two initially independent 
condensates begin to overlap. For £ (t) being nonzero, we see clearly from the coefficient b d that the interference term 
in (J7J gives a contribution to the density expectation value. Based on Eqs. (fTTjl. (fT2")) and (fT3|) . we have 



b d (i) \C{t)\{Nx-i) 



c d (i) 



4+1 



(15) 
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FIG. 2: Based on the numerical calculations of Eqs. J5), (|9} and (|10[) . shown in Fig. 2(a) is the relation between bd/cd and (, 
while shown in Fig. 2(b) is the relation between bd/cd and Ni = N2 for £ = 0.001. It is shown clearly that the interference 
term can play an important role in the density expectation value for |£| being larger than iVf 1 . 



and 

mo = \gt)\(N 2 -z) ^ 

ad (i) i + 1 

It is easy to understand that when Ni |£ (t)\ > 1 and N2 |C (^)l > 1, &d can not be omitted, and thus there would be 
clear interference patterns. This shows clearly that for large particle number, a small |£ (t)\ can play important role 
in the density expectation value when |£ (t)\ > iV-f 1 and |C(*)I > ^2^ 1 - 

Generally speaking, even in the presence of the interatomic interaction, |£ (t)\ is much smaller than 1 because </>i</>2 
is an oscillation function about the space coordinate. As shown above, however, a nonzero value of £ (t) can give 
significant contribution to the density expectation value for large iVi and N%. In Fig. 2(a), we give the relation 
between bd/cd and £ for TVi = iV2 = 10 3 based on Eqs. ^ and (jTU|) . The relation between bd/cd and iVi = N2 for 
C = 0.001 is shown in Fig. 2(b). Generally speaking, for iVi |C| >> 1 and iV 2 |C| >> 1, one has bd/^/adCd ~ 1. 
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III. THE EVOLUTION EQUATIONS OF THE SYSTEM 

For the many-body wave function ^^^2) in the presence of the interatomic interaction, the evolution equation can 
be obtained based on the standard quantum mechanical principle. After the double-well potential is removed, the 
evolution equation is given by 

in — Q- — = H f ^ NlN2 , (17) 
where the Hamiltonian of the whole system in the first-quantization method is 

N ti 2 N 

= E - 2^ V ? + 9 E - rj) ■ (18) 

i—l i<j 

Here the coupling constant g — ATrh 2 a s /m with a s being the scattering length. In the above expression of the 
Hamiltonian, we have used the well-known two-body pseudopotentials. 

It is well known that the action principle is quite useful to derive the GP equation for a single condensate (See 
for example Ref. Q). Similarly, we consider here the evolution of '^n 1 n 2 based on the action principle. To get the 
evolution equation based on the action principle, we first give the general expression of the overall energy of the whole 
system. After removing the double-well potential, the overall energy of the whole system is 



E = j y* NlN2 H f y NlN2 d i r 1 ---d ^n. (19) 

After straightforward derivations, the exact expression of the overall energy is given by 

E = E km + E int , (20) 

where the kinetic energy Ekm is 

f ( N ti 2 \ 

Ekin = J y* NlN2 [Y, ~^ V H y^N^r! ■ ■ ■ d 3 r N 

+^__ e -*e V0 . . V ^ + £^! V ^ . V A . (21) 
Ira im ) 

In addition, the interaction energy Ei nt of the whole system is given by 

JV 



E int = J V* Ni n 2 ( ffE S ( r * ~ r ^ I ^N 1 N 2 d 3 r 1 ■ ■ ■ d 3 r N 



I IdV 



+Rc (h 4 |^| 2 ftfae** + h 5 (0I) 2 4> 2 2 e 2i ^ + h 6 |0 2 | 2 <f>*ihe 1 ^)] , (22) 
where the coefficients are given by 

1 ^ i!t!(JVi-*-2)!(iVa-*)!' 

mH NuN 2 - 2) Al Nl \N 2 \\tF 

h 2 - 2^ 



(Nt - i)\ (N 2 - i - 2)V 

i—Q 

min(Afi-l,JV 2 -l) , 2i 

h 3 = 4.1-Ai!A 2 !|c| 



iHl{N 1 -i-iy.{N 2 -i-iy. 
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ft 4 = 2^ 



i\ (i + 1)1 (Ni — i — 2)1 (N 2 — « — 1)! ' 

2^!^! |C| 2<+2 



4=0 



i\ (i + 2)\ (JVi - i - 2)! (N 2 -i- 2)1 ' 
4AlN 1 \N 2 \\C\' 2 



6 ^ i!(i + l)!(JVi-i-l)!(JV 2 -i-2)!" 1 J 

By using the ordinary action principle and the energy of the whole system, one can get the following coupled 
evolution equations for fa and fa: 

^ = (24) 
dt NxSfa' v ' 

lH -dt = NlW, (25) 
where SE/Sfa and SE/Sfa are functional derivatives. 

IV. THE CROSSOVER FROM |C| « A^ 1 TO \(\ » N' 1 
A. the case of |f| << TV" 1 

In Sec. II, we have shown that the nonzero interference term in the density expectation value origins from the 
assumption that C (t) can be a nonzero value after the overlapping between two initially independent condensates. 
In the last section, we have given the evolution equations about fa and fa. Thus, an important question emerges 
naturally: is it physical for £ (t) being nonzero with the development of time? 

Based on the evolution equations (f!M)) and (|23|) given in the last section, one can understand easily that ( (t) 
becomes nonzero after the overlapping between two initially independent condensates in the presence of the interatomic 
interaction. Before removing the double- well potential, iVx |Cj = and AT 2 |£| = because there is no overlapping 
between fa and fa. After removing the double- well potential and at the beginning of the overlapping between two 
condensates, Ni |£| << 1 and N 2 |£| << !• For N\ \(\ << 1 and N 2 |C| << lj the overall energy can be approximated 
as 

ti 2 r 

E w - — / dV (JViV^I • Vfa + N 2 Wfa ■ Wfa) + 
2m J 

| J dV^N 1 (N 1 -l)\faf + N 2 (N 2 ~l)\fa\ 4 + AN 1 N 2 \fa\ 2 \fa\ 2 y (26) 

The above overall energy is obtained by setting £ = in Eq. ([20|) . 

Based on Eqs. (|24|) and (|25|) . the approximate evolution equations for fa (r, t) and fa (r,t) are then 



^dfa h 2 9 , , 

^ + (27) 



and 



where 



and 



iK^=-^V 2 fa + V 2 fa, (28) 
at 2m 



V^i^-^glfal 2 + 2N 2 g\fa\ 2 , (29) 



V 2 = (N 2 -l)g\fa\ 2 + 2N ig \fa\ 2 



(30) 
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We see that V\ is not equal to V 2 for g 7^ 0. This leads to an important result that £ (t) can be nonzero when there 
is an overlapping between two initially independent condensates. One can understand this result further through the 
following equation which determines the evolution of £ (t) : 

ifi^ = J f{v,t)MW, (31) 

where / (r, t) is a nonzero function given by 

/ (r, t) = (N 2 + 1) g |0 2 1 2 - (N, + 1) g |^| 2 . (32) 

The above analyses show clearly that why £ (t) becomes nonzero after the overlapping between two initially inde- 
pendent condensates for j ^ 0. Obviously, for g being zero, £ (t) — at any further time because d( (t) /dt = 0. 
Our numerical calculations in the following section also show that C (t) can be a nonzero value in the presence of the 
interatomic interaction. 

In the presence of the interatomic interaction, it is inconsistent to assume that C (t) is always zero with the devel- 
opment of time. If ( (t) = 0, the evolution equations Eqs. (f2~4"|) and ([23]) are exact. However, based on these two 
evolution equations, |£| will increase from zero after the overlapping. 



B. the case of \(\ » N' 1 

As shown above, the evolution of 4>i and (f> 2 is determined by the expression of the overall energy and the evolution 
equations (f2"4"|) and (|2"5|) . It seems that there is no exact solution for these evolution equations considering the fact that 
even there is no exact solution for three-dimensional nonlinear Schrodinger equation. However, for |£| >> N^ 1 and 
ICI >> N% , we find that there is a quite simple evolution equation by introducing an effective parameter order. This 
would also contribute to our understanding of the interaction-induced coherence process for two initially independent 
condensates. 

For the cases of N± |C| >> 1, -ZV2 jC! >> 1 and ^1 ~ ^2, first we introduce the following effective order parameter 
<i> e (r, t) which is given by 

$ e (r, t) = v/A^>i (r, i) + yfihj*' <f>2 (r, t) . (33) 

Based on this effective order parameter, the density expectation value can be approximated well as 

n d {v,t) ~ **(r,t)$ e (r,t) 

= a'J0i(r,i)| 2 + 26^xRe[ e ^^(r,i)0 2 (r,i)] +c^|0 2 (r,t)| 2 , (34) 

where 

a' d = N x , 

b' d = ^Njh, 

c' d = N 2 . (35) 

To compare with the exact expression of the density expectation value given by Eq. Fig. 3 shows the ratio 
Ai = dd/a'd, ^2 = b d /b d and A 3 = Cd/c' d . We see that for Ni\(\ >> 1, N 2 \(\ >> 1, the approximate density 
expectation value (r, t) <!> e (r, t) agrees very well with the exact expression of the density expectation value. 

Based on the effective order parameter, for JVi |C| >> 1, N 2 |C| >> 1 and N\ ~ N 2 , we also find that the overall 
energy of the whole system can be approximated very well as 

E'a PP — E' kin + E int , (36) 

where 

E'kin = 4" / v$ : ■ v$ edV : (37) 
2m J 

and 



#L = f / dV\<f e \ 4 



(38) 
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FIG. 3: Shown is the values of Ai, A2 and A3. We see that for A''! |£| >> 1, N2 |C| >> 1 an d Ni ~ N2, the density expectation 
value given by the effective order parameter agrees well with the exact expression of the density expectation value given by Eq. 



For the cases of Ni \(\ » 1, 7V 2 | CI >> 1 an d Ni ~ N2, it is easy to verify that E' kin Ek% n based on the analogous 
analyses about the density expectation value. For Ni |£| >> 1, N 2 |C| >> 1 and N% ~ N%, one can also prove the 
result of E' int w E int . Based on Eq. (f3"5)) . E' int can be expanded as: 



where 



TP' 



9 
2 

+Re 



dV 



+ /3 2 |<? 


i 2 | 4 + /?3^l| 2 |02| 2 




+ ft W) 2 <&f?* c + ft |0 2 | 2 #6>e** 


A = 


Nl 


ft = 


Nl 


03 = 




ft = 




ft = 


2N 1 N 2 , 


ft = 


4JV 2V /iViJV 2 . 



(39) 



(40) 
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FIG. 4: Shown is the relation between /ij//3j (i = 1, 
A^2 |C| >> 1, hi w 1, which means that .Ej nt « I?mt. 



,6) and £ for Ni = N 2 = 10 3 . We see that for ATi |<| >> 1 and 



To compare with the exact expression of the interaction energy given by Eq. ([22]) . Fig. 4 shows the relation between 
hi/fa (i = 1, • • • ,6) and C for N r = N 2 = 10 3 . It is shown clearly that for N% |C| >> 1 and N 2 \(\ >> 1, h/fa w 1, 
and thus ^' nt ps S int . 

In this situation, the overall energy can be approximated well as 



EL 



2m 



V$! • \7<5> e dV 



9 I dV\$<* 



(41) 



The evolution of the effective order parameter can be then obtained by using the action principle. Based on this 
approximate energy, after removing the double-well potential, it is quite interesting to note that the evolution of the 
effective order parameter can be described very well by the ordinary Gross-Pitaevskii equation: 



ih- 



dt 



2m 



V 2 $ e +.g|$ e | 2 $ e . 



(42) 



Based on Eqs. and (|4"2"|) , we see that the emergence of the effective order parameter $ e gives us strong evidence 
that the coherence is formed in the interaction process between two initially independent condensates, and thus results 
in the emergence of high-contrast interference fringes. The effective order parameter and the approximate density 
expectation value suggest strongly that a full coherence is formed between two initially independent condensates for 
the cases of Ni |C| >> 1 and N 2 |C| >> 1- 
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V. THE EVOLUTION OF THE DENSITY EXPECTATION VALUE ACCORDING TO THE 

EXPERIMENTAL PARAMETERS 

Now we turn to give the theoretical results of the density expectation value according to the experimental parameters 
in Ref. I5J where clear interference patterns were observed for two initially independent condensates. In the experiment 
of Ref. 5(, N = 5 x 10 6 condensed sodium atoms were confined in a magnetic trap with u x — 2tt x 18 Hz, u> y = uj z = 
2ir x 320 Hz. A blue-detuned laser beam of wavelength 514 nm was focused into a light sheet with a cross section 
of 12 /zm by 67 /im. The long axis of the laser beam was perpendicular to the long a;— axis of the condensate. For a 
laser power of 14 mW, the barrier height is about 1.4 /iK, which is much larger than the chemical potential fj, = 0.03 
fiK. For this laser beam, the two condensates can be regarded to be independent because they are well separated and 
the tunneling effect can be omitted. With these experimental parameters and the s— wave scattering length a s — 2.75 
nm, the initial profile of the two condensates is shown in Fig. 5. At the initial time, the overlapping between two 
condensates can be omitted safely, and thus £ (t = 0) = 0. After the double-well potential is removed, the evolution 
of the density expectation value rid-x (x, t) — J rid ( r > t) dydz (in unit of N/2) is given in Fig. 5 through the numerical 
calculations of Eqs. (f24|) . ([25|) and We see that there is a clear interference pattern in the density expectation 
value which agrees with the experimental result. Shown in the inset is the evolution of |£| for these parameters. We 
see that |£| increases from zero with the development of time. For the expansion time of 40 ms, the numerical result 
of |C| shows that bd/cd ~ 1- To check further that the nonzero value of \C\ does not originate from numerical error, 
with the same initial conditions, we have verified in the numerical calculations that |£| is always zero (smaller than 
10 -10 ) if the scattering length is assumed as zero. 

Based on the evolution of the density expectation value shown in Fig. 5, we see that the overall width (about 300 
/im) of the ultracold gases in the x— direction does not increase obviously. This is due to the fact that the initial 
density distribution is cigar-shaped, and thus the expansion in the x— direction is very slow, while the expansion in y 
and z directions is much quick. In the experiment of Ref. 19], one can see clearly that there is no obvious expansion 
in the long x— axis for cigar-shaped condensate. For two initially independent condensates, when the double-well 
potential is switched off, one should note that in the regime close to x = 250 /im shown in Fig. 5, the ultracold 
gases expand rapidly in the a;— direction because in this regime the ultracold gases have higher kinetic energy. Thus, 
although the total width of the system in the x— direction does not increase obviously, the rapid expansion in the 
central regime leads to the overlapping between two initially independent condensates, and results in the interference 
effect. The overall width of the system shown in Fig. 5 is smaller than the experimental result of about 500 /jm [f|. 
This difference may come from the expansion of thermal cloud in this experiment [f|. After 40 ms expansion, the 
numerical result in [l(| for two coherently separated condensates also showed that the overall width of the system is 
about 300 /Ltm. 

VI. THE EVOLUTION OF THE DENSITY EXPECTATION VALUE FOR DIFFERENT COUPLING 

CONSTANTS 



We see that the interatomic interaction plays an essential role in the emergence of the interference effect for two 
initially independent condensates. Generally speaking, increasing the particle number will enhance the effect of the 
interference term in the density expectation value. Based on Eqs. (j2"4")) and (|25p. increasing the coupling constant g 
has the effect of increasing £ (t). Together with the relation between bd/cd and ( illustrated in Fig. 2(a), this shows 
that increasing the interatomic interaction will enhance the effect of the interference term. To show more clearly the 
interaction- induced coherence process between two initially independent condensates, in this section we consider the 
density expectation value for different coupling constants. 

We consider here the evolution of the density expectation value for one-dimensional case. At t = 0, to give 
a general comparison, the initial wave functions for two independent condensates are assumed to be identical for 
different coupling constants. The initial wave functions are respectively given by 



01 (X[,t 



0) 




1 



exp 



(xi - Xii) 
2A2 



(43) 



02 (Xl,t 



0) 




1 



exp 



(xi - xu) 
2A| 



(44) 



In the above wave functions, wc have introduced a dimensionless variable xi = x/l with I being a length. In the 
present work, we assume that Ai = A2 = 0.5 and 2^2 — %ii = 4.5. For these parameters, at t — 0, the two condensates 
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0.003 




Position (jiim) 



FIG. 5: Shown is the evolution of the density expectation value calculated using the experimental parameters in Ref. [f|. It is 
shown that there is a clear interference pattern after the overlapping between two independent interacting condensates. Shown 
in the inset is the evolution of |£ (t)\ for the experimental parameters in Ref. [f|. 

are well separated. In the numerical calculations of the coupled equations given by Eqs. (|24"|) and (f2"5| . it is useful to 
introduce the dimensionless variable r = Erf/H with Ei — ti 2 /2ml 2 , and dimensionless coupling constant <?; = N\g/Eil. 
In addition, the particle number is assumed as N\ = N2 = 1.0 x 10 5 . In real experiments, interatomic interaction 
plays very important role in the initial ground-state wave function of the condensates. However, in principle, one can 
prepare the state given by Eqs. (|43p and (144)) by adjusting the trapping potential for different coupling constants. In 
this section, the identical initial wave functions for different coupling constants would be helpful in the comparison of 
the density expectation value for different coupling constants. 

For t > 0, we consider the evolution of the density expectation value in free space. The evolution of <pi and <p2 is 
obtained based on the numerical calculations of Eqs. (f!M|) and ([^5]) . From tfii and <j>2, we can get £, and thus the 
density expectation value based on Eq. ([7]). Shown in Fig. 6 is the evolution of the density expectation value for 
different coupling constants. It is shown clearly that increasing the coupling constant has the effect of enhancing the 
coherence effect, and results in higher contrast in the interference patterns. 
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FIG. 6: After two initially independent condensates are allowed to expand freely, shown is the evolution of the density 
expectation value nd(xi,r) (in unit of Ni + N%) for different coupling constants gi. Shown in the inset of each figure is 
the relation between bd/cd and dimensionless time r. For two ideal condensates shown in Fig. 6(a), we see that there is no 
interference pattern even there is an overlapping between two condensates. For the case of gi — 1 shown in Fig. 6(b), we see 
that low-contrast interference patterns begin to emerge due to the interaction-induced coherence process. In Fig. 6(c) and Fig. 
6(d), we see that there are high-contrast interference patterns. In particular, in Fig. 6(d) for gi — 20, two initially independent 
condensates can be regarded to be fully coherent near r = 1. 



VII. THE CALCULATIONS BASED ON THE SECOND-QUANTIZATION METHOD 



The essential reason for the emergence of the interference term of two initially independent condensates lies in 
that because of the exchange symmetry of identical bosons and interatomic interaction, the two initially independent 
condensates become coherent after the overlapping between the two condensates. The physical mechanism of this 
interaction-induced coherence can be understood further based on the second quantization method. Thus, in this 
section, we give the density expectation value and evolution equations based on the second-quantization method. The 
merit of the second-quantization method lies in that the indistinguishability of identical bosons is satisfied when the 
correct commutation relation of the field operator is used. 
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A. the density expectation value 



For two initially independent condensates comprising particle number N± and N 2 , the corresponding quantum state 
is (See for example Refs. 

\N lt N a ) = -^&r&)»>\0), (45) 

where 5„ is a normalization constant to assure (iVj., iVjj | iVi , N 2 ) = 1. a\ (a 2 ) is a creation operator which creates a 
particle described by the single-particle state fa (fa) in the left (right) condensate. Similarly to Eq. the quantum 
depletion is omitted in this sort of quantum state. 

For two initially independent condensates, for the state given by Eq. (|45p . it seems that there is no interference 
term in the density expectation value with the following simple calculation: 

n d (r,t) = (N 1 ,N 2 ,t\^\N u N 2 ,t) 

= (N 1; N 2 ,t\ (a+Si |0! (r, i)| 2 + «|S 2 | (r,t)| 2 ) |JVi,JV 2 ,t> 

+2 x Re ([Ni,N 2 ,t\ a\a 2 \N u N 2 ,t) fa (r, t) fa (r, i)) (46) 
= N x \fa(v,t)\ 2 + N 2 \fa(v,t)\^ . (47) 



In the above equation, the field operator ^ is expanded as "J = a\fa + a 2 fa + • • • with <Si and a 2 being the annihilation 
operators. One should note that, to get (|47|) from (146)) 
This holds when / fa (r, t) fa 2 (r, t) dV = 0. When 

(the last term in Eq. (|4"6")l ) is zero in (r, t). 

As shown in the previous sections, fa and fa will become non-orthogonal in the presence of the interatomic 
interaction. The operators a\ and a 2 can be written as 



there is an implicit assumption that a\ and a\ are commutative. 

0, it is easy to understand that the interference term 



ai,a 2 



a ± = I VfadV, (48) 

and 

a 2 = I ^fadV. (49) 



Here ^ is the field operator. By using the commutation relations of the field operators [ty (ri, i) , \1/ (r2, t)\ = and 
[ty (ri, t) , ^ (r 2 , t)] = 5 (1*1 — r 2 ), it is easy to get the commutation relation 

[Si,aJ]=C*. (50) 

We see that <5i and a 2 are not commutative any more for J fafadV being a nonzero value. In this situation, it is 
obvious that one can not get the result (|4"T|) from (|4"6"|) any more. This means that one should be very careful to get 
the correct density expectation value for two initially independent condensates. 

It is well-known that the field operator should be expanded in terms of a complete and orthogonal basis set. 
Generally speaking, the field operator ^ can be expanded as: 

$ = aifa + k<f>' 2 -\ , (51) 

where fa and fa are orthogonal normalization wave functions. Assuming that fa = (3 (fa +afa) 7 based on the 

\ j 2 

conditions J fafa 2 dV = and / \fa\ 2 dV = 1, we have \/3\ = (l - |C| 2 ) and a = - (*■ Based on 

% = J *(^)*dV, (52) 

we have 

% 

a 2 = — + Cai. (53) 
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It is easy to get the following commutation relations: 



[k,k] = fe f ] = 0, \k,%i] = 1, 
[01,01] = [a\,a[] = 0, [Si, a}] = 1, 
[fc,ai] = [fc,a}]=0. 



(54) 



Because k and a{ are commutative, it is convenient to calculate the density expectation value rid (r, t) by using the 
operators k and a\. After straightforward derivations, the exact expression of the density expectation value is 



rid 0, t) 



iVi!iV 2 ! v 

+ k<t>' 2 



k _ 



N 2 



(8!) 



2Vi 



fct 



N 2 



+ Ca\ 



|0) 



where the coefficients arc 



a d |0! (r, i)| 2 + 2/3 d x Re (e*°fi (r, i) 2 (r, t)) + 7d \fa (r, i)| 2 
»2iV 2 !(7V 1+l -l)W 1 (l-|C| 2 ) A ' 2 "|C| 21 



/3d = 



Id 



E 

i=0 
JV 2 - 

E 

i=0 
AT 2 - 

E 

i=0 



m 


(iV x - 1)! (JV a - i)! 






/ ,\JV 3 -i- 

M)!(l-|C| 2 ) 


_1 IC| 2l+1 


il (i + 1 


)! (iVx - 1)! (iV 2 - * 


-1)! ' 


iS2jV 2 !(Ar x - 


/ \N 2 -i- 

M)!(l-|C| 2 ) 





ililNxl (N 2 - i - 1)\ 



In addition, the normalization constant s„ is determined by 

/ 



iV 2 !(7V 1 +i)!(l- |C| 2 ) ^ ' id 24 



i=0 



i!i!JVi! (JV 2 — i) 



= 1. 



(55) 

(56) 
(57) 
(58) 



(59) 



The above density expectation value is obtained based on the second quantization method. Although it seems that 
the expressions of the coefficients given by Eqs. ([56]) , (f57|) . (|58|) and (|59|) are quite different from the results calculated 
from the many-body wave function, we have proven that ad = ad, [3d — bd, Jd — Cd and S„ = A n . Thus, the density 
expectation value given by Eq. (|55p is the same as the result calculated from the many-body wave function \& n 1 n 2 
which satisfies the exchange symmetry of identical bosons. 



B. the evolution equations 

Based on the second quantization method, after removing the double- well potential, the Hamiltonian in the second- 
quantization method is 



H s = [ dv( — V$ f • V$ + ^$t$t§$ 
J \2m 2 

After straightforward calculations, the overall energy of the whole system is 

E = Ekin + Ei n t , 

where the kinetic energy E^m is given by 



(60) 



(61) 
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E ki n = J (N u N 2 ,t\^V& ■^\N 1 ,N 2 ,t)dV 



= J d V 



2m 



2m 



• V0i + ^-v$ • V</) 2 



2m z 2m 

In addition, the interaction energy i?j nt of the whole system is given by 

= | "1 l^ll 4 + "2 |02| 4 + «3 |0l| 2 |02| 2 

+Re (a 4 |0i | 2 bXfae*' + a 5 (<^) 2 <fie 2i *<> + a 6 
where the coefficients are given by 



ai 



a 2 



a 3 



a4 



a 5 



a 6 



JV 2 

E 

j=0 



5^iV 2 !(JV 1 +t-2)!iVi(iV 1 -1) 
i!t! (JVi - 2)! (JV 2 - *)! 



(i-ICI 2 ) 



Nn—i 



ici 2 \ 



llp c 



N f? 5*N 2 \(N 1+ i)\ 
^ i\i\N l \(N 2 -i-2) 

i=0 
JV 2 -1 



o - if 2 ) 



N 2 -i-2 



ici 2 \ 



4~? n N 2 \{N 1 +i- ly.Ni 



(iVi - 1)! (N 2 - i - 1) 
4S2 7V 2 ! + i _ i)!iv : 



E 

i=0 
iV 2 -l 

E 

i=0 
N 2 -2 

E 

i=0 
N 2 -2 

E »!(i + i)!(iVi-l)!(JV2-t-2) 



JVo-i-l 



icr 



t!(i + l)!(iVi-2)!(JV2- 


t - 1)! 


2Z 2 n N 2 \(N 1+ i)\ 




i\ (i + 2)! (Ni - 2)! (JV 2 - 


t-2)! 


4E 2 n N 2 \ (JVi + 1)! 





^wHci*) 



ATo-i-1 



ICI 



2i+l 



N 2 -j-2 



ICI 



2i+2 



]V 2 



ICI 



2i+l 



(62) 



(63) 



(64) 



We have proven that the above overall energy is equal to the results based on the many-body wave function, by 
checking that hi = a* for i = 1, 2, • • ■, 6. 

When the non-orthogonal property between <f>\ and cf) 2 are considered, because the derivations of the density 
expectation value and overall energy are quite cumbersome, the same results based on the first-quantization method 
and second-quantization method give us strong evidence that our derivations are correct. 



VIII. QUANTUM FLUCTUATIONS AND ORTHOGONALITY OF THE WHOLE QUANTUM STATE 

In both the first-quantization and second-quantization methods of the previous calculations, the quantum depletion 
originating from the elementary excitations at zero temperature is omitted in the quantum state of the whole system. 
Based on the Bogoliubov theory of the elementary excitations, for the initial quantum state, the number of particles due 

3 /2 

to the quantum depletion is of the order of (a/I) and thus the quantum depletion is negligible for Bose condensate 
in dilute gases considered in the present work. With the development of time, the role of quantum depletion can also 
be omitted. Because of the factor e ±ik r (|k| is the wave number of the elementary excitations) in the wave function 
0k of the elementary excitations, a simple analysis shows qualitatively that (0i|0k) and (02|0k) are of the order of 
|£| e -(l k l£) with L being the spatial size of the system. This exponential decay of (0i|0k) and (0 2 |0k) originates from 
the integral where there is spatially oscillating phase factor in the wave functions of the elementary excitations and 
condensates. Thus, the contribution to the effective order parameter and density expectation value due to elementary 
excitations can be omitted safely. 

In the preceding paragraph, it is shown that in calculating the evolution of the density expectation value, the role 
of quantum depletion can be omitted. However, the quantum depletion plays important role in the consistency of our 
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theory. Here we discuss mainly the orthogonality for the quantum state of the whole system, especially about the 
subtle problem that whether the nonorthogonality between fa and fa violates the orthogonality of the whole system 
which must be satisfied. 

This sort of problem about orthogonality exists also for a single condensate. For a single condensate, if the quantum 
depiction is omitted, the many-body wave function is 

V(r 1 ,---,r N ,t) = 4>(r 1 ,t)---4>(T N ,t). (65) 

Based on the many-body Schrodinger equation and the action principle, it is easy to get the following GP equation 

= "£ V20 + Vext (r ' t ) ( t ) + 9{N-l) H 2 4>- (66) 

Assume that there are two orthogonally initial quantum states ^ a (i = 0) = fa\ (ri,t = 0) • • • <f>A (*n, t — 0) and 
*b (t = 0) = ((>b (ri, t = 0) • • • <j>g ( r N,t = 0). With the development of time, because of the nonlinearity of the GP 
equation, <f>A and 4>b may become non-orthogonal in the presence of interatomic interaction. It is obvious that 
/ ^* A ^Bdri ■ ■ ■ dr N = [J (j>* A (r, t) cj> B (r, t) dr] . Although for large N, this integral can be approximated as zero, 
based on the consideration of consistency, the physical mechanism of the exact orthogonality between Si! a and ^> B is 
an interesting problem. 

We think that the effect of quantum depletion omitted in Eq. (|63|) is the essential reason for this sort of unphysical 
nonorthogonality. For the Hamiltonian (fP8| of the whole system, the term g Y^zj $ ( r * — r j ) ^ s non-factorable about 
the coordinate, and thus the exact solution of the many-body wave function is non-factorable too. For a single 
condensate, the exact many-body wave function can be assumed as 

* 'exact (ri, ■■■,r N ,t) = aj>(ri,t) ■ ■■<j)(r N ,t) F N (n, ■■■,r N ,t). (67) 

Here Fn accounts for the non- factorable component, and c is a normalization constant. In the hard-sphere approx- 
imation, Fjy = if | r i — r j | < 2r for any i ^ j (r is the hard-sphere radius), and F/v = 1 for otherwise situation. 
Fjv represents the quantum depletion, and was successfully used to calculate the quantum depletion of superfluid 
liquid 4 Hc in Ref. [20(. The omission of the quantum depletion means that we approximate Fn as 1. For dilute Bose 
condensed gases, r « I. Thus, the wave function can describe well the condensate. When the orthogonality 
of the whole system is considered, we stress here that the wave function ^ exact should be used. The non-factorable 
factor Fn assures the orthogonality of the many-body quantum state, because VP exact is the exact solution of the 
many-body Schrodinger equation. For dilute Bose condensed gases, these analyses lead to two results: when the 
dynamic evolution of dilute Bose condensed gases is considered, omitting the quantum depletion can give us quite 
good description; when the consistency especially the orthogonality is considered, one should consider the role of 
quantum depletion. 

It is natural to generalize the above analyses to two initially independent condensates. The nonorthogonality 
between <j>i and 4>2 does not mean in any sense the violation of the orthogonality of the quantum state of the whole 
system. We stress here again that when orthogonality is considered, we should check whether the quantum state 
of the whole system satisfies the orthogonality. When the quantum depletion is considered, the many-body wave 
function takes the following form 

^WiNa (ri, r 2 , ■■■ ,r N ,t) = c^ NlN2 (r 1 ,r 2 , ■ ■ ■ ,r N ,t)F N (r x , ■■■ ,r N ,t) . (68) 

Here ^>n!N 2 is given by Eq. ([5]). ^n 1 n 2 satisfies exactly the many-body Schrodinger equation, and thus 
this sort of wave function satisfies the orthogonality condition. Assume that is constructed by <pi 

and cf>2, while ^ e N 1 N 2 ^ s constructed by cf>f and <p B . Assume further that 0^, <p A , <fif and <p B are orthog- 
onal to each other at the initial time, so that v E'^^v 2 and ^^^ 2 are orthogonal at the initial time. Omit- 
ting the quantum depletion and roughly speaking, the integral between ^^^r 2 and ^^^ 2 is of the order of 
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f((tf)*<f>fdT J (4>t)* 0fdr /(</>£) *</>fdr /(^)*</>fdr (here a + (3 + 7 + 6 = N). We see that similarly 

to a single condensate, the orthogonality between and 4'^"^, is quite good for large N when the quantum 

depletion is omitted. But to assure the exact orthogonality of the whole quantum state, we must consider the role of 
quantum depletion. Although we do not provide here a fully rigorous proof, these analyses lead to the result that: 
although for dilute Bose condensed gases, the omission of the quantum depletion will not play important contribution 
to the evolution of (f>i and (f>2, it's the quantum depletion that assure the orthogonality of the whole quantum state. 

Up to date, we mainly use the bases <p\ and <p2 to study the dynamic evolution and density expectation value. Con- 
sidering the nonorthogonality between fa and fa in the presence of interparticle interaction, it is attractive to discuss 
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the physical picture using the orthogonal bases 4>i and (j)' 2 introduced in Sec. VII. If we use the orthogonal bases (pi and 
<t>' 2 , the quantum state of the whole system is \N U N a ) ~ (a\) Nl 0/f3+Cal) N2 |0) ~ Em=o C m (a[) Nl+m (k^) N2 - m |0). 
We see that the number of particles in the orthogonal modes <f>\ and <j> 2 are no more definite. This quantum state 
becomes a superposition of different number of particles in the orthogonal modes <f>i and (f>' 2 . This result is very 
natural because the coupling (interaction) between two initially independent condensates leads to the coherent trans- 
fer of particles between two condensates. In Ref . [2l[ , it was shown that for a single condensate, a coherent state 
description of the Bose condensed system is a robust state in the presence of the interactions between the condensate 
and its environment. In the present work, our research shows that the interaction between two condensates makes 
each condensate become a coherent superposition of different particle number. The results in Ref. [2l[ about the 
robustness of the coherent state imply that £m=o C m (a\) Nl+m (k^) Na - m |0) is a robust quantum state even when the 
coupling with the environment exists. 

IX. SUMMARY AND DISCUSSION 

In summary, upon expansion, we calculate the density expectation value of two initially independent condensates. 
It is found that there is a nonzero interference term when the interatomic interaction and the exchange symmetry 
of identical bosons are both considered carefully. In fact, it is well-known that the interaction plays an essential role 
in the formation of the order parameter of Bose-condensed gases, i.e. the formation of a stable coherent property. 
Here, we provide an example in which the interaction induces coherent evolution between two initially independent 
condensates. 

In the present popular viewpoint, the high-order correlation function and quantum measurement process are used 
to interpret the interference patterns observed for two initially independent condensates. Although we find that there 
is already high-contrast interference patterns in the density expectation value for the experimental parameters in Ref. 
Q, it is still possible that the high-order correlation function and thus the mechanism of the interference due to the 
measurement process play an important role in the emergence of the interference patterns. Thus, to investigate more 
clearly the interference pattern due to the measurement process alone, we believe an experimental investigation of two 
independent ideal condensates would be very interesting, because there is no interference term for ideal condensates 
in the density expectation value. In the last few years, the rapid experimental advances of Feshbach resonance where 
the scattering length can be tuned from positive to negative make this sort of experiment be feasible. 

Note added: The present paper is an expansion of our two previous unpublished works [22l. |23|. In addition, after 
the submission of the present paper, we noticed a paper [2~ij by L. S. Cederbaum et al., where the role of interaction 
in the emergence of interference patterns for two initially independent Bose condensates is stressed. 
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